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INTRODUCTIOri 


Considerable  progress  has  been  made  toward  a  comprehensive 
understanding  of  the  formation,  development  and  microphysical 
structure  of  clouds  and  fogs  formed  over  the  summit  of  Great 
Dim  Fell  together  with  their  related  absorption  of  radiation 
at  optical  and  infra-red  wavelengths. 

The  experimental  program  incorporates  measurements  of 
condensation  nuclei  concentrations  and  associated  meteorological 
parameters  made  within  the  air  feeding  into  these  clouds, 
utilizing  a  mobile  laboratory,  in  addition  to  microphysical  and 
transmission  measurements  made  within  such  clouds  at  the  summit 
of  Great  Dun  Fell.  Equipment  available  for  these  studies  which 
has  been  developed  or  modified  for  ground-based  usage  includes 
Kelly  probe  electrostatic  disdrometers  and  two  PMS  spectrometer 
probes  (FSSP)  for  the  measurement  of  cloud  droplet  distributions 
together  with  a  modified  Mee  CCN  counter  and  a  PMS  ASASP  probe 
for  measuring  the  smaller  aerosol  particles.  A  passive  optical 
scattering  instrumen-fc  of  novel  design  for  the  simultaneous 
measurement  of  visibility  and  cloud  liquid  water  content  has 
been  developed  and  utilized  in  field  studies  at  Great  Dun  Fell. 
Associated  measurements  of  the  transmission  of  infra-red  radia¬ 
tion  at  about  11um  have  been  made  with  a  3arnes  transmissometer. 

A  comparison  of  measured  extinction  coefficient  values 
for  radiation  of  this  wavelength  with  values  calculated  from 
microphysical  measurements  is  provided  in  Section  III. 

Theoretical  work  has  included  the  development  of  an 
airflow  model  capable  of  predicting  the  microphysical  properties 
of  clouds  over  the  summit  of  Great  Dun  Fell  from  a  knowledge  of 
the  relevant  parameters  concerning  the  air  parcels  from  which  the 
clouds  are  formed.  In  addition,  computational  studies  have  been 
made  of  the  role  of  entrainment  of  undersatura ted  air  upon  the 
droplet  size  distribution  for  situations  relevant  to  the  field 
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investigations.  This  investigation  of  turbulent  mixing  is 
described  in  Section  I  of  this  report.  Section  II  provides 
a  detailed  comparison  of  the  theoretical  cloud  model  with 
field  observations  made  in  differing  atmospheric  situations. 

The  work  presented  in  this  report,  ujhilst  still  in  its 
preliminary  stages  and  capable  of  much  further  refinement,  has 
nevertheless  shoun  that  the  microphysical  properties  of  clouds 
formed  over  Great  Dun  Fell  may  be  predicted  from  a  knowledge  of 
the  parameters  of  the-  air  belouj  cloud  base  with  a  notable 
degree  of  accuracy.  Furthermore,  the  initial  investigations 
of  the  transmission  of  infra-red  radiation  through  these  clouds 
suggest  that  this  parameter  may  be  directly  related  to  their 
microphysical  properties.  Thus  the  prospect  of  being  able  to 
estimate  the  attenuation  characteristics  likely  to  arise  during 
the  formation  and  development  of  clouds  and  fogs  from  the 
preceding  meteorological  conditions  appears  increasingly 
realistic. 


I.  TURBULENT  MIXING  MODEL 
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Ue  have  made  some  calculations  to  investigate  the  influ¬ 
ence  of  turbulent  entrainment  of  environmental  air  upon  the 
adiabatic  cloud  droplet  spectrum.  Lie  assume  the  turbulent 

mixing  to  be  roughly  represented  by  diffusion,  uiith  a  constant 

2  - 1 

coefficient  K  (m  s  )  determined  on  dimensional  grounds  from 
the  original  spectral  scale  X  (m)  of  the  entrained  parcel,  and 
€  (ms  ),  the  rate  of  dissipation  of  turbulent  energy  uithin 
the  cloud.  It  should  be  noted  that  this  approximation  is  best 
suited  to  the  case  uihere  \  is  much  greater  than  the  character- 
istic  eddy  size  yithin  the  clouds:  K  =  (X^e)^.  Neglecting 
buoyancy  effects  and  assuming  the  turbulent  motions  transport 
sensible  heat,  liquid  water  and  water  vapour  simultaneously, 
we  have 


^(R,t) 

fit 

-  K7^T(R,t)  =  Tgyap(R.t) 

(1) 

-  Kv2L(R,t)  *  L.g^gp(R,t) 

(2) 

-  Kv^^^CR.t)  =  P3^3p(R,t) 

(3) 

and 


■-evap^^-^^  =  -‘*^j’^n(r,R,t)^ 

[S(R,t)-^  +  |3]dr 

(4) 

(5) 

« 

a  (R,t)  =  -L  (R,t) 

*^evap  — ’  evap  — ’ 

(6) 

(7) 

j)g(R,t)  =  1.a.1Q~^exp(~|-^gP-^) 

(g  m"^) 

(3) 

n(r,^,t)dr  is  the  number  (m~^)  of  drops  at  (R_,t)  with  radii  in 
the  interval  r,r+dr,  and  a  is  a  diffusion  length  we  have  put 
equal  to  5jim.  .A=  ZM^^^/RT^y,  B  =  3vM^mg/4njp^Mg ,  r  is  the  droplet 
radius,  and  the  molecular  weights  of  water  and  the  aerosol. 


V  the  number  of  ions  per  aerosol  molecule,  q  the  surface 
tension  of  the  aqueous  solution  drop  against  air,  the 
density  of  water,  L  the  latent  heat  of  evaporation  and  R  the 
universal  gas  constant.  Ue  assume  that  all  droplets  have 
identical  nuclear  mass  m^.  Uhile  these  equations  could  be 
made  more  precise  by  improving  the  expressions  for  J3,pg  and 
the  droplet  growth  equations,  such  changes  were  deemed 
unnecessary  at  this  stage  in  view  of  the  very  approximate  relation¬ 
ship  between  the  diffusion  problem  and  the  mixing  process  we 
are  studying. 

The  numerical  computations  assume  that  a  spherical  blob 
of  radius  \  is  placed  at  t=Q  at  the  centre  of  a  spherical  cloud 
shell  of  outer  radius  Because  buoyancy  can  be  neglected 

spherical  symmetry  is  preserved  throughout  the  diffusion  phase 
change  process.  Ue  assume  the  boundary  conditions 

5i-(l,t)  =  ^(l,t)  =  ilv(l,t)  =  f^(l,r.,t)  =  0 
6t  5l  6l  ^  ^ 

at  1  = 

T(0,t),  L(0,t),  j3y(0,t)  and  n(Q,rj,t)  are  finite,  for  all  r^  and  t. 

In  the  calculation,  diffusion  and  phase  change  proceed  sequentially. 
That  is,  at  the  beginning  of  each  time  step  all  parameters  are 
known  at  each  spatial  grid  point  Ij^.  The  grid  spacing  chosen 
was  Al=  ^=\/10,  where  lj^  =  (i=1)ul.  Diffusion  proceeds  via 

a  forward  differencing  algorithm  with  the  right-hand-sides  of 

n 

the  first  three  equations  equal  to  zero  for  a  time  ut=CAl)  /6M 
set  for  numerical  stability.  Then  all  motion  stops  and  phase 
change  (evaporation  or  condensation)  takes  place  at  each  grid 
point  according  to  equations  (£*)  to  (3).  This  continues  for  a 
time  ilt  after  which  the  time  variable  is  advanced  by  ut  and  the 
diffusion  step  repeated.  The  time  step  for  the  phase  change 
equations  was  set  equal  to  O.IAt  and  was  always  less  than  O.ls. 


The  stabilitv  of  the  diffusion  algorithm  was  checked  for 
each  run  by  comparing  the  results  of  the  diffusion  calculation 
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in  the  absence  of  phase  change  with  the  analytical  solution 
for  the  spherically  symmetric  problem  with  the  zero  flux  boundary 
condition.  For  all  runs  considered  the  numerical  and  analytical 
solution  agreed  to  within  Q.15i 
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ABSTRACT 

A  description  is  giv/en  of  experiments  performed  to 
investigate  the  effect  of  the  entrainment  of  undersaturated  air 
upon  the  microphysical  structure  of  cap  clouds  at  Great  Oun 
Fell,  Cumbria.  The  results  obtained  are  compared  with  the 
predictions  of  simple  theoretical  models.  Three  case  studies 
are  presented;  in  case  I  the  properties  of  the  cloud  uere 
essentially  adiabatic;  II  provided  evidence  for  mixing- inter¬ 
mediate  betu/een  the  homogeneous  and  extreme  inhomogeneous 
types;  in  III  the  mixing  was  extremely  inhomogeneous.  Con¬ 
sistent  ujith  this  picture,  the  time  constant  for  turbulent  erosion 
of  the  entrained  parcels  of  air  is  significantly  greater  in 
case  III  than  in  case  II.  Spectral  changes  predicted  by  a 
model  of  turbulent  entrainment  (Baker  and  Latham,  1981)  agree 
well  with  observation. 


I.  INTRODUCTION 


Baker  and  Latham  (1979)  and  Baker,  Corbin  and  Latham 
(1980)  have  presented  calculations  which  suggest  that  the  influ¬ 
ence  of  dry  air  entrainment  on  the  microphysical  structure 
of  a  cloud  is  strongly  dependent  on  the  ratio  of  the  time 
constants  6r  turbulent  diffusion  of  entrained  air,  and 

for  the  evaporation  of  drops  in  response  to  the  under sa turation 
T^.  Measurements  of  microphysical  properties  of  clouds 
enveloping  Great  Dun  Fell  (GDF)  in  Cumbria  described  by  Blyth 
et  al  (1980)  suggested  that  entrainment  has  an  appreciable 
effect  on  these  clouds.  In  addition, the  indications  were  that 
for  the  clouds  sampled  the  effect  of  the  mixing  is  closer 
to  that  of  the  extreme  inhomogeneous  description  of  Baker 
and  Latham,  )  t  than  the  classical  description, 

which  implies  that  «  "I* 

In  this  paper  three  detailed  case  studies  are  discussed 
which  throw  additional  light  on  the  effect  of  entrainment 
upon  the  properties  of  clouds  over  Great  Dun  Fell. 

In  case  I  entrainment  effects  appeared  to  be  negligible  and 
cloud  was  thus  essentially  adiabatic.  In  cases  II  and  III 
the  liquid  water  contents  were  substantially  sub-adiabatic, 
indicating  that  significant  entrainment  of  undersaturated 
environmental  air  occurred.  In  each'  case  a  simple  dynamical  and 
microphysical  model  of  the  cloud  is  used  to  relate  the  micro- 
physical  measurements  to  .the  meteorologicel  conditions.  An 
attempt  is  made  to  assess  the  influence  of  entrainment  upon 
the  cloud  properties  in  terms  of  a  simple  diffusive  model  of  the 
mixing  process  developed  by  Baker  and  Latrn-m  (1951). 
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2.  INSTRUMEIVTATIQIM  AND  EXPERIMENTAL  PROCEDURES 

The  basic  arrangements  and  procedures  uiere  similar  to 
those  employed  by  Blyth  et  al  (1980)  so  uill  not  be  reproduced 
here.  Lie  merely  outline  the  instrumentation  and 
draui  attention  to  additional  devices  employed. 

At  the  summit  of  the  hill  drop  spectral  data  uiere  obtained 
using  the  Keily  Probe  (HP)  (Corbin  et  al  1978),  cases  I  and  III, 
and  a  PMS  Axially  Scattering  Spectrometer  Probe  (ASSP)  (Ryder 
1976),  case  II,  and  for  a  small  period  of  case  III.  Occasional 
measurements  of  drop  concentration  and  size  distribution  uiere 
made  using  a  pulsed  Laser  Holography  System  (Conuiay  et  al  '1980,1951 
A  high  frequency  acoustic  sounder  (Moulsley  and  Cole  1979) 
uias  employed  to  monitor  the  height  of  the  lowest  inversion 
base,  taken  as  cloud  top  height  (see  e.g.  Caughey  et  al  1978, 

1981). 

The  ASSP-100  is  regarded  as  primarily  a  droplet  sizing 
instrument.  A  comparison  was  carried  out  between  droplet 
spectra  from  the  holographic  laser  system  and  the  ASSP-1D0 
(see  Conway  et  al  1981).  Good  agreement  was  obtained  for  both 
mean  radius  and  spectral  shape.  Droplet  concentrations  (and 
hence  liquid  water  dontent)  are,  however,  known  to  be  over¬ 
estimated  by  our  ASSP.  Comparison  with  droplet  concentrations 
deduced  from  the  holographic  technique  during  the  present 
studies  resulted  in  a  correction  factor  of  0.5  for  the  ASSP 
droplet  counts  and  this  has  been  applied  to  the  whole  data 
set.  This  figure  is  in  good  agreement  with  that  deduced  from 


other  studies  of  ASSP  measured  liquid  water  contents  and  those 
from  other  instruments  (Kitchen,  1979,  Slingo  et  al  1981). 


The  results  obtained  uiith  the  Keily  probe  and  the  corrected  PM£ 
ASSP  data  uere  compared  in  detail  during  the  period  1635-1650 
on  15  May  (Case  III).  The  cloud  properties  shoued  no  o\/erall 
trend.  Figure  1  shouis  a  comparison  oF  droplet  spectra  and  a 
comparison  of  parameters  is  shouin  in  table  1.  The  error  terms 
represent  the  short  term  variability  in  the  data.  It  is  seen 
that  the  agreement  is  good. 

Observations  of  cloud  condensation  nucleus  spectra  (CCiM) 
were  made  Just  below  cloud  base  using  a  modified  Mee  CC^'i 
device.  Occasional  profiles  of  wind  speed  and  direction  at 
2m  above  the  local  terrain,  and  of  temperature  were  taken 
from  just  below  cloud  base  to  the  mountain  top. 


ID 

3.  ADIABATIC  MODEL  OF  CLOUD  OROPLEf  EVyOLUriDIM 

The  model  described  in  this  section  calculates  the  cloud 
droplet  size  distribution  oF  an  air  parcel  subjected  to 
adiabatic  ascent.  Initially  the  air  has  a  relative  humidity 
of  39%  with  a  known  size  distribution  of  droplets,  calculated 
from  the  CCN  activity  spectrum,  in  vapour  equilibrium.  Ue 
start  with  a  humidity  of  99%,  sinc'e  for  higher  humidities 
the  large  nuclei  do  not  have  enough  time  to  attain  their 
equilibrium  size. 

The  equations  for  the  model  are  Very  similar  to  those 
used  by  Lee  and  Pruppacher  (1977)  without  entrainment,  the 
only  significant  difference  being  that  the  droplet  growth 
equation  has  been  adapted  so  that  no  specific  knowledge  of 
aerosol  types  or  mass  is  required.  The  equations  were  inte¬ 
grated  numerically  using  a  simple  forward  difference  method. 

For  droplet  growth  we  have 

S  -  ^  ^ 

_  r  r-' 


r .  dr 


dt 

p  RT 
w 

L^w 

iT 

e(T)0(r)M,, 

L  W 

K(r)T 

Ur  -  y] 

(1) 


where  A  =  2M^o-/RTp^,  B  =  ,  S  is  the  super sa t ur a t ion  , 

r  the  droplet  radius,  and  the  molecular  weights  of  water 

and  the  aerosol,  the  number  of  ions  per  aerosol  molecule, 

O’  the  surface  tension  of  the  aqueous  solution  drop  against  air, 

the  density  of  water,  e(T)  the  saturation  vapour  pressure 

with  respect  to  a  flat  water  surface  at  temperature  T.  D(r)  and 
respectively 

K(r)  are| dif f usivi ty  of  water  vapour  and  thermal  conductivity 
of  air  at  a  drop  radius  r,  L  the  latent  heat  of  evaporation 
and  R  the  universal  gas  constant. 

The  time  variations  of  the  vapour  mixing  ratio, w^,  and 
liquid  water  content,  w^^, 
sedin^ntatlnn)  by 


are  given,  for  adiabatic  ascent  (no 
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du^ 

JT 

du^ 

dT 


du^ 


(2) 


(3) 


where  n  is  ills  density  of  air  and  N(r)  the  concentration  of 
^  8 

droDs  of  radius  r.  The  temperature  variation  is  taken  to  be 
-dT  ^  ^  ^  L  dw^ 


dt 


dt 


dS 

P_ 

dw 

V  - 

-  (1  +  S) 

>L 

f— ) 

+ 

dt 

0  e 

dt~ 

[RgT^ 

Ut  J 

R  T 
a 

li)  is  the  updraught  speed,  Cp  the  specific  heat  of  dry  air  at 
constant  pressure  and  g  is  the  acceleration  due  to  gravity. 
The  time  variation  of  supersaturat ion  is 

„  1 

(5) 

where  P  is  the  pressure,  P  =  i^s  ratio  of  the  molecular 

weights  of  dry  air  and  water,  Rg  is  the  specific  gas  constant 
of  dry  air. 

Finally  the  pressure  variation  is  given  by 

jPU! 


dP 

dt 


R  T 
a 


(6) 


We  now  show  how  input  data  for  the  model  were  prepared 
from  the  CC(\I  activity  spectrum  and  obtain  the  required  form 

of  the  droplet  growth  equation.  The  equations  of  Kohler 

« 

(Pruppacher  and  Klett  1978)  relate  the  equilibrium  size  of 
a  drop  with  nucleus  mass  m^  to  the  ambient  supersaturation  S. 


r 


(7) 


A  and  B  are  as  previously  defined.  Differentiation  of  this 
equation  to  find  a  maximum  in  5  yields  a  critical  radius: 


^C  \A  / 

R-;'jriting  (7)  usin.i  (B)  we  have 


(8) 


A 

r 


1  r 


S 


(9) 
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The  Mee  activity  spectrum  gives  the  number  of  drops  activated 

at  a  given  supersaturation .  We  divided  the  total  concentration 

of  drops  into  categories  whose  size  varies  from  about  UOcm”^ 

to  Q.05cm“^  as  the  drop  activation  supersaturation  (critical 

supersaturation)  decreased.  This  ensured  a  small  range  of 

critical  supersaturations  for  all  drops  in  each  category,  and 

enabloj  us  to  assign  one  value  for  the  critical  supersaturation 

(S  ) ,  to  all  drops  in  one  category  with  only  small  error.  It 

was  subsequently  found  that  changing  the  size  of  the  categories 

of  the  CCN  caused  no  significant  change* in  the  calculated  adiabatic 

droplet  spectrum.  For  the  1^^  category  the  critical  radius 

(r  is  related  to  the  assigned  critical  supersaturation  (S  ). 
c  i  c  i 

by  setting  r  .=  (r  ).  in  equation  (9)  giving 

c  * 


(r  )  =  l-±- 

^^c^i  3  (S^) 


(10) 


C  1 


The 

radius  of  droplets  at  cloud  base  is  obtained  by  substituting  the 
cloud  base  supersaturation  in  (9).  We  assumed  a  supersaturation 
S  =  -0.01  giving  r^^  as  a  solution  of  the  equation 


,0.01  =  ^  -  1  ^^c>i  ^ 

^i  3 


(11) 


The  input  data  can  now  be  obtained;for  category  i  we  have  a  drop 
concentration  AN.  with  critical  radius  (r  ),  from  (10),  and  cloud 
base  radius  r^  from  equation  (11). 

Finally  we  rewrite  the  droplet  growth  equation  using  (8) 

For  the  category  i: 


ri  ^ 

dt 


S  _  ^  i  <^0^' 

r.  3  - 3- 

i  ri-’ 


(12) 


)0(r)M..  M(r)T  \RT  "  V 


e(T)0( 


w 


We  have  eliminated  8  which  depends  on  m  and  can  use  our  knowledge 


of  fna  critical  radius 


Mast  atmospheric  aerosol  particles  are  mixed;  i.e.  they 

consist  of  uiater  soluble  and  insoluble  parts,  (ile  have  shoun 

in  the  appendix,  that  the  above  equations  can  still  be  used 

for  these  mixed  particles  with  only  very  little  error. 

speed 

The  updraught/liiCj: )  of  horizontal  distance  x  from  the 
summit  uas  taken  to  be 

Iil(x)  =  f  grad(Jf)  UCje)  *  (13) 

where  U('*)  is  the  local  wind  at  2  metres  above  the  ground  -  averaged 
over  about  5  minutes  -  and 

^stimated  from  measurements  made  at  different  locations  up 
the  hill  using  our  instrumented  Land  Rover,  grad(Jf)  is  the 
local  gradient  estimated  from  a  map  scale  (1:25, OQQ),  and  f 
is  a  factor  introduced  to  allow  for  the  retarding  effect  of 
friction  at  the  ground.  Assuming  an  approximately  logarithmic 
wind  profile  we  took  f  =  1.5.  Subsequent  analysis  using  the 
linear  theory  for  turbulent  flow  over  a  hill  (Jackson  and  Hunt 
1975)  suggests  this  to  be  a  reasonable  estimate.  Values  of 
temperature  and  pressure  required  for  the  calculations  were 
measured  at  cloud  base. 


4.  DROPLET  LOSS  TO  GROUND 


The  effect  of  loss  of  drops  to  ground  by  eddy  diffusion  and 
gravitational  sedimentation  uas  investigated  using  the  equation 


/,,  /  v^NCrK  .  y  iM. 

Vt  ■  az  bz  ^  az 


5l\l(r) , 


aN(r) 


(14) 


uhere  N(r)  is  the  drop  concentration  of  radius  r,  is  the 
terminal  velocity  and  K(z)  is  the  eddy  diffusivity  at  height  z. 
The  equation  was  solved  numerically  using  a  .forward  difference 
method  and  starting  with  a  typical  drop  size  distribution  which 
was  the  same  throughout  the  cloud.  The  integration  was  carried 
out  for  the  time  it  takes  a  parcel  of  air  to  travel  from  the 
cloud  base  to  the  mountain  top. 

The  primary  difficulty  in  carrying  out  these  calculations 
is  in  the  choice  of  the  eddy  diffusivity,  K.  We  used 


K(z)  =  □.4zu«  z  <  15Qm 
h(z)  =  60  u*  z  >  150m 


(15) 


where  u,  is  the  friction  velocity  deduced  from  2m  and  10m  winds 
at  mountain  top.  The  first  expression  for  K  is  based  on  the 
assumption  of  near  neutral  stability,  and  the  second  simulates 
a  highly  diffusive  outer  layer.  This  is  certainly  a  large  over¬ 
estimate  for  any  situation  in  which  mechanical  turbulence  is 
dominant  (eg.  Bradley  1900).  Convection  and  the  occurrence  of 
rotary  type  motions  may  both  act  to  increase  deposition;  however  no 
information  is  available  on  the  magnitude  of  these  effects.  It  is 
felt  that  the  above  eddy  diffusivity  will  still  lead  to  an  upper  boun 
for  droplet  loss  to  ground*.  For  the  upper  boundary  condition 
K  was  assumed  to  fall  rapidly  to  zero  at  cloud  top,  i.e.  at  the 


base  of 


the  capping  inversion.  Support  for  this  notion 


comes  from  turbulence  structure  observations  in  the  vicinity  of 
inversions  over  level  country  (e.g.  Caughey  and  Palmer  1979) 
and  from  the  appearance  of  an  echo  free  region,  implying  small 
turbulence  levels,  above  the  inversion  in  the  acoustic  records 


5 


The  results  showed  that  depletion  was  only  significant  for 
clouds  of  vertical  depth  <  20Qm.  The  cloud  depth  is  a  minimum 
above  the  mountain- top.  We  estimate  the  depth  to  be  greater  than 
5QQm,  between  50  and  IQQm,  and  around  3Q0m  for  the  cases  I,  II 
and  III  respectively.  Hence  absolute  values  of  drop  concentration 
and  liquid  water  content  must  be  treated  with  some  caution, 
especially  in  case  II,  where  subadiabatic  liquid  water  contents 
may  not  be  due  to  entrainment  alone.  The  preferential  deposition 
of  larger  drops  by  gravitational  sedimentation  was  negligible  for 
a  typical  range  of  drop  sizes,  and  caused  no  significant  change 
in  the  drop  size  distribution  in  the  time  required  for  the  cloud 
to  flow  from  cloud  base  to  the  mountain  top. 
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5.  THE  FIELD  EXPERIMENTS 
Case  I  10  May  1979 

A  deep  iaolated  cap  cloud  uias  located  over  the  ridge  at 
GDF  and  precipitation  uas  falling  into  it  from  an  initially 
much  higher  deck  of  altostratus  associated  uith  an  approaching 

_  A 

warm  front.  Winds  uiere  south  easterly,  about  10m  s  at  10m 

above  the  summit.  A  section  of  the  acoustic  record  for  the 

period  ~11.3Q  to  12.10  BST  is  given  in  Figure  (2).  This 

illustrates  the  presence  of  fairly  weak  thermal  plumes  extending 

up  to  ~50-10Qm  in  height  from  the  hill  top.  It  is  also  apparent, 

from  the  absence  of  an  elevated  layer  echo,  that  no  significant 

entrainment  interface  existed  within  the  lowest  3D0n,  and  this 

tends  to  support  the  vieu'  that  on  this  occasion  en  t  r  a  i  nme''!  t 

effects  at  the  mountain-top  were  likely  to  be  quite  small.  (The 

thin  stationary  echo  at  ~90m  is  from  one  of  the  communication 

towers  on  GOF ,  )  Figure  3a  shows  that  after  the  formation  period  the- 

adiabe  tic 

cloud  had  an  approximately  adiabatic  liquid  water  content.  (The^ 


value  was  calculated  from  the  height  and  temperature  of  cloud 
base  using  the  model.)  Figure  3b  shows  the  varieition  in  15 
minute  averages  of  the  drop  concentration  and  P  the  percentage 
of  3-secand  drop  counts  which  satisfy  one  of  the  conditions 

N  >  N  -i-  2  /N 

□  r  N  <  N  -  2  /¥  ^ 

as  a  measure  of  the  small  scale  inhomooeneity  in  the  cloud,  as 

(1990). 

used  by  Blyth  et  al^-^  W  is  the  number  of  drops  in  the  3-S9cond 

histogram  and  W  is  the  15  minute  average  drop  concentration. 

Care  was  taken  to  establish  that  the  number  of  histograms 

satisfying  one  of  the  criteria  (16)  was  not  s i gn i f i ce n tl y 

affected  by  trends  in  the  mean  drop  concentration  over  a  15 

minute  period.  The  Iom  veli.'Es  of  P  indicate  a  fairly  ho-cr-necu?  do. 
e  t  r  If  c  t  u  r  e  . 


Figure  4  shows  a  typical  spectrum  in  the  cloud,  and  the 
adiabatic  model  spectrum  calculated  for  the  same  time. 

Figure  3  and  U  strongly  suggest  a  cloud  with  near  adiebatic 
properties,  whose  nearly  homogeneous  structure  indicates  little 
effect  from  interaction  with  its  environment.  Examination  of 
anemograph  records  from  the  mountain  top  shows  that  the  air 
was  at  least  as  turbulent  in  case  study  I  as  in  cases  II  and  III. 
It  is  therefore  reassuring  for  the  interpretation  of  cases  II 
and  III  that  P  is  small,  indicating  that  water  loss  to  ground  is 
not  introducing  structure  into  the  cloud  on  the  scales  examined. 

_  <1 

The  rainfall  rate  in  the  area  of  GDF  did  not  exceed  2mm  hr 
and  rough  calculations  confirm  that  the  effect  of  scavenging 
of  droplets  by  precipitation,  or  the  breakup  of  rain  drops  had 
a  negligible  effect  on  the  droolet  spectra. 


Case  II  Ik  May  1979 

The  cap  cloud  formed  in  a  weakly  convective  boundary  layer 

capped  by  a  strong  subsidence  inversion;  routine  ascents  showed 

that  this  was  generally  at  around  600m  above  lowland  areas,  but 

it  was  forced  to  rise  in  the  vicinity  of  the  hill,  and  was  observed 

by  the  acoustic  sounder  to  be  between  lOm  and  IDOm  above  the 

mountain  top.  Cloud  base  was  60Cm  below  the  mountain  top  at 

07.00  aST  and  thereafter  rose  steadily  to  be  150  below  the  summit 

at  15.00  0ST.  The  10m  wind  at  the  summit  was  from  210°  with  a 
_  1 

speed  of  ~10m  s  .  Figure  (5)  presents  a  section  of  the  acoustic 
recording  between  1030  and  1120  0ST.  The  very  dark  elevated 
layer  echo  which  varies  in  height  from  ~1Cm  to  above  100m  is 
from  the  subsidence  inversion.  Undulations  were  present  on  a 
range  of  scales,  i.e.  from  a  few  hundred  nietres  to  mesoscele 
V'  ;  i  c- 1  i  ori  r>  of  "10-1:. ''.n  (e.g.  the  din  which  occurs  t;etiu-'n  10. 


ia 

and  11.00  BST.  The  dark  region  extending  from  the  inversion 
into  the  near  surface  region  probably  represents  entrained  air 
mixing  into  the  boundary  layer.  Thus,  uhilst  entrainment  was 
almost  certainly  occurring  locally  the  intensity  of  the  echo 
from  the  inversion  (reflecting  the  inversion  strength)  suggests 
that  it  would  be  weak.  The  sounder  also  indicates  the  presence  of 
weak  convection  under  the  inversion.  7he  CCIM  distrib¬ 

ution  was  much  more  maritime  than  in  either  case  I  or  case  III. 
Curve  P  of  Figure  6  shows  the  percentage  of  3-secand  drop 
counts  which  satisfy  the  criterion  (  10.  Although  the  mean 
liquid  water  content  was  substantially  subad iaba t i c ,  especially 
early  in  the  period, (see  Figure  7,  which  shows  15  minute  averages 
of  liquid  water  content),^  the  regions  of  highest  drop  counts  had 
liquid  water  contents  which  approached  the  adiabatic  value 
estimated  from  the  observations  ofcloud  base  height,  and  drop 
concentrations  similar  to  those  predicted  by  the  model.  However 
the  measured  concentrations  must  be  treated  with  caution  in 
view  of  the  adjustment  described  earlier.  Further  the  droplet 
spectrum  observed  for  these  regions  agrees  quite  well  with  the 
calculated  adiabatic  spectrum,  especially  early  in  the  day 
as  can  be  seen  from  Figure  8.  The  spectra  in  the  high  and  low 
drop  concentration  regions  (defined  by  (  16))  are  shown  in 
Figure  9.  The  shapes  of  the  spectra,  taken  with  the  proximity 
of  the  dry  air  source,  suggest  that  the  inhomogene i t  ies  in  the 
cloud  were  produced  by  the  entrainment  of  undersaturated  air, 
and  this  is  supported  by  the  acoustic  sounder  trace,  (Figure  5). 

The  turbulent  mixing  model  outlined  by  Beker  end  Lathem 
(1931)  we  s  utilised  to  predict  changes  to  the  high  con  l  e  ri  t  r  a  t  i  on 
spec  true-.  (A)  in  Fi-,;  ure  9  r>"ultirm  f  roei  the  mixing-in  of 
env  i ;  T,  a  ]  air.  A  sh'c'.fir  was  to  sf  u'hc  tn.;  the  Ip'.' 

Orj^*'_.’'ltJr''inri  -iLii  (c)  U]‘*0  The  .-"'w’.i  O' 


parameters  utilised  in  these  calculations,  based  on  the 
meteorolooical  observations,  uiere:  e,  (the  turbulent  enercy  , 

3  9  p  3 

dissipation  rate),  2  x  TO  to  2  x  10  ms  ,  estimated 

from  the  mean  wind-speed  at  an  altitude  of  10m,  the  assumption 

of  neutral  stability  and  an  assumed  surface  roughness  of  2cm; 

and  X  (the  radius  of  the  spherical  entrained  blob)  =  50  to  IDOm, 

estimated  from  the  acoustic  sounder  traces.  Ue  assumed  that 

the  su pe r sa t ura t ion  S,  estimated  from  ascent  data,  was  around 
- 1 5 

-0.2,  m^  =  10  g,  and  that  the  temperatures  before  mixing  of 
the  cloudy  and  environmental  air,  and  T^,  were  50  (measured) 
and  30  (estimated)  respectively.  The  spectrum  A  (L  =  D.^45g  m~^) 
was  subdivided  into  5  size-categories. 

Ourve  0  in  Figure  9  is  calculated  on  the  basis  of  this 
mixing  model,  and  is  seen  to  be  in  close  agreement  with  the 
observed  spectrum  3,  which  has  a  reduced  liquid  water  content 
of  Q.OSg  m"^.  As  the  mixing  process  proceeds,  in  the  calcu¬ 
lations,  there  are  many  regions,  at  particular  times,  where 
L  =  0.8g  m“^,  but  examination  snows  that  the  droplet  spectra 
are  always  similar  to  0.  It  was  found  also  that"  this  spectrum 
was  insensitive  to  reasonable  variations  in  the  input  parameters 
(Baker  and  Lathem,  1931).  This  agreement  between  observed  and 
calculated  spectra  encourages  us  to  believe  that  the  mixing 
model  is  capable  of  describing,  with  adequate  precision,  the 
influence  of  entrainment  upon  the  droplet  spectra. 

Examination  of  curve  P,  Figure  6,  shows  that  there  is 
an  entrainment  maximum  around  11.00  after  which  a  continual 
decrease  occurs  until  near  the  end  of  the  period.  Despite  this 
the  low  count/high  count  spectra  showed  tfiat  the  form  of  the 
mixing  re'nained  as  described  above.  The  i  nt-iomagene  i  ty  at  the 
erid  of  tiie  p‘  i  iod  war  di:e  tc  patcby  clrvic!  ae  bt-.rv  ror. 

above  the  mcnrit-^ln  too. 


The  large  rise  in  drop  concentrations  towards  the  end  of 

adiabatic 

the  period  was  partially  reproduced  by  thej^model.  The  increase 
in  the  calculated  value  of  N  was  due  to  higher  updraught  speeds 
near  cloud  base  as  this  moved  up  the  hill.  This  produced 
higher  peak  supersaturations  and  so  activate^  more  nuclei. 

However  the  narrowing  of  the  spectrum  was  greater  than  predicted. 
This  may  be  a  consequence  of  the  crudity  of  the  model  but  would 
result  from  a  change  to  a  more  continental  CCN.  Measurements  with 
the  Mee  device  were  unavailable  to  test  this  latter  possibility, 
but  a  change  to  a  more  continental  CCM  had  occurred  by  the  next 
day . 


Case  III  15  May  1979 

The  day  was  characterised  by  light  SU  winds,  (less  than 
tm  s”  at  the  mountain  top).  A  general  cover  of  stratus  and 
stratocumulus  cloud  enveloped  the  mountain.  At  11.00  BST  cloud 
base  was  500m  below  mountain  top  and  rose  to  first  clear  the 
mountain  at  1A00  BST.  There  was  a  weak  subsidence  inversion 
(much  weaker  than  in  case  II)  about  150m  above  the  mountain  top. 
This  was  clearly  revealed  in  the  acoustic  sounder  record,  a 
section  of  which  is  given  in  Figure  (10).  The  important  features 
to  note  are  that  the  echo  layer  is  much  weaker  than  on  the  l^ith 
and  is  also  highly  contorted  on  a  range  of  scales.  Convection 
beneath  the  interface  is  also  much  stronger  than  on  the  previous 
day.  These  features  suggest  a  si-uation  in  which  entrainment 
effects  are  likely  to  be  prominent.  Regions  of  entraining  air 
can  be  seen  extending  from  tiie  interface  to  near  the  ground  at 
-Id.lO  aST  and  again  et  --1/.  .20  BST. 

During  the  ptriod  u:i  lo  It*. 30  the  adiabatic  liqjid  water 
cont-.r,t  fell  B.y  ;i;  ■  : -j f  (Figure  1lX  eltnciugh  initial]'; 
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!:■  jnsport  af  moisti.'rp  p.”d  CC-'j  in  the  L.Tjn::-ry  layer.  After  IP.'jn 
tile  acoustic  sounder  s>''Oueri  po  t  r  3  in;ii.en  t  into  thp  cloud  throuijh 
tha  cajipinQ  in\/prsion  on  scales  of  seueral  tens  O'  mpt.es,  (Ficure 
IQ).  A  larqe  fall  in  mean  drop  count  occurred  as  the  cloud 
became  very  inhomogeneous,  as  shoun  in  Figure  (12),  but  the  peak 
drop  concentrations  observed  in  the  patches  of  thickest  cloud  re¬ 
mained  unchanged  until  about  13.10.  During  this  time  the  high 
and  louj  count  spectra  (derived  by  the  criteria  16)  became  very 
similar  in  shape  (Figure  13A,3)  suggesting  that  the  mixing  uas 
extremely  inhomogeneous.  Uithin  this  period  the  observed  speci-rum 
OBcame  appreciably  broader,  (Figure  1A,A).  Houever,  it  was  no 
broader  than  the  model  adiabatic  spectrum  calculated  for  the 
samp  time  using  CCD  measurements  made  at  11. AO.  Curve  C  in 
Figure  13  is  the  spectrum  predicted  from  the  turbulent  mixing  mcdel 
for  input  parameters  estimated  from  the  mateorolocical  measurements 
in  the  same  manner  as  For  Figure  9.  Their  values  were:  e  = 

5  X  IQ'^m^s"^,  \  =  50m,  S  =  -0.2,  =  9°C,  =  s'^C.  The 

w  C 

initial  spectrum  A,  (L  =  0.2Ag  m  was  suodivided  into  12  size- 
categories.  There  is  seen  to  be  good  agreement  between  the  observed 
and  calculated  spectra,  9  and  C  respectively. 


During  the  remainder  of  the  experiment  cloud  base  remained 
near  the  mountain  top.  The  curve  of  the  number  of  histoorams 
outside  the  2  PJ  criterion  shows  a  marked  lull  in  mixing  around 
16.10  3ST.  This  was  confirmed  by  the  acoustic  sounder.  This 
period  was  characterised  by  a  narrower  spectrum  as  shown  in 
Figure  1A(a).  These  experiments  therefore  provide  further  strong 
evidence  for  the  importance  of  understanding  the  nature  of  the 
entrainment  processes  in  any  study  of  the  microphysical  development 
of  real  clouds. 
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6.  DISCUSSION 

The  general  conclusions  from  these  three  case  studies  may 
be  interpreted  in  terms  of  the  recently  developed  ideas  on  the 
interaction  of  clouds  with  their  environment  (Baker  et  al, 

1930).  For  extreme  inhomogeneous  mixing  the  time  constants 

for  eddy  dissipation  and  droplet  evaporation,  my/T£  ®,  whereas 
for  classical  mixing  "*  0- 

In  case  I  local  turbulence  levels  were  high,  but  the  cloud 
was  deep  and  exhibited  an  approximately  adiabatic  structure. 

Any  entrainment  effects  produced  small  changes  at  the  hilltop 
sampling  point.  Ue  may  estimate  the  relative  extents  to  which 
the  entrainment/ mixing  processes  in  cases  II  and  III  are  inhomo¬ 
geneous  by  considering  the  ratio  ^  t  j/t^  )  j  j  j/ C  -j-/ - ^  )  j  j  .  Ignoring 
the  contribution  to  the  turbulence  kinetic  energy  budget  from 
buoyancy  production  we  assume  that  most  of  the  turbulence 

energy  is  generated  through  mechanical  production.  The  above 
ratio  can  then  be  written  es  follows: 


(t 

^III 
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m  i  X  i  n  0 

lengths,  a  (''5um)  the  scale  length  associated  with  the  conden¬ 
sation  coefficient  and  S  the  undersaturation  of  the  entrained 


air.  The  acoustic  sounder  traces  indicate  that 

the  droplet  data  show  that  rj^~12um,  rj^j~Fj;n,  and  we  employ  the 
values  of  e  and  S  pres  on  ted  earlier:  to  40, 

t  there  fore  folio;.'’;  that 


I 


Buoyancy  production,  as  indicated  by  the  acoustic  sounder,  was 
certainly  much  larger  during  case  III  than  case  II  and  uiould 
therefore  hav/e  acted  to  reduce  the  above  ratio.  Houever,  inclusion 
of  this  factor  would  not  he  expected  to  alter  the  conclusion  that 
the  above  ratio  is  consistent  with  the  experimental  observations 
that  the  cloud  is  much  more  inhomogeneous  in  case  III  than  in  case 
II. 

A  further  prediction  of  the  inhomogeneous  mixing  model  of 
Baker  et  al  (1930)  is  that  the  drop  concentrations  of 

the  largest  drops  are  greater  than  those  obtained  by  purely 
adiabatic  growth.  Comparisons  of  calculated  adiabatic  spectra  and 
observed  mean  spectra  show  no  clear  evidence  for  this  (the  in¬ 
accuracies  in  the  calculated  updraught  profiles  do  not  seriously 
affect  this  result  since  the  measured  CCN  distribution  leads  to 
calculated  spectra  much  more  similar  in  shape  to  the  broadest 
rather  than  the  narrowest  observed  spectra  for  a  range  of  likely 
updraughts).  In  case  II  the  lack  of  broadening  can  be  explained  if 

entrainment  effects  are  confined  close  to  the  region  near  the  mountain 
top  where  the 

^wind  shear  across  the  inversion  interface  is  greatest;  droplets 
unaffected  by  the  entrained  air  will  not  now  have  time  to  grow  to 
sizes  which  are  significantly  in  excess  of  those  resulting  from 
adiabatic  growth.  In  case  III  we  observed tha t  the  highest  con¬ 
centrations  of  large  d ro ps  oc cn r red  in  regions  where  the  degree 
of  inhomogeneity  and  trie  liquid  water  content  were  highest. 
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THE  EFFECT  OF  MIXED  NUCLEI  ON  THE  DROPLET  GROUITH  EQUATION 

It  is  assumed  that  the  uater  insoluble  part  is  a  sphere 
of  radius  x.  The  Kohler  equation  for  the  equilibrium  drop 
size  is  (Pruppacher  and  Klett,  1978) 

B 


S 

T 


-  ;c3 


(A1) 


where  A  and  B  are  as  defined  in  the  body  of  the  paper, 


Differentiating  as  before  gives 


B  = 


(r^3  _ 

3r  ^ 


(A2) 


and  rewriting  (A1) 


n  /  3  3.2 

S  =  ^  ^  ^ 

r  - 


( r  ■ 


(A3) 


3r , 


Setting  r  =  r^^  the  critical  radius  is  related  to  the  critical 


supersaturation  by 


P  .  2  ^(1  ^ 

^  3  3 


(A4) 


This  equation  is  identical  to  (ID)  for  completely  soluble 
aerosol  if  ix^/2v^^)  «  1 

Pruppacher  and  Klett(1973),  table  6.3,  present  some  values 

related  to  this  parameter.  They  use  salt  and  silicon  dioxide 

respectively  as  the  soluble  and  insoluble  parts.  With  these 

compounds  the  error  is  less  than  3'4  for  5  and 

1  s 

m  >  10“  c.  where  m  is  the  mass  of  the  insoluble  part  and 

TJ  -  ’  X 

II 

m.  the  total  nuclear  mass.  Work  by  Hanel  (1976)  has  shown 
N 

that  all  CCN  likely  to  be  activated  at  Great  Dun  Fell  satisfy 
thesi'  conditions  well. 

For  ttif  equilibriun  droo  radius  r  at  cloud  base  we  have 


-D.D1  =  - 


-  :c^)2 


3. 

(r  -  X  } 


3r. 
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If  uje  can  assume  (— )  t  (— )  «  1,  uje  can  use  first  order 

c 


approximations  to  obtain 
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3r‘ 


:c  3  ^3 

1  _  2(— )  ^  (£) 
t-  c  ^  . 


(A7) 


This  equation  is  identical  to  (11)  if  ue  let  (r  ) 

c-'mixed 


^^c^soluhle  introduces  only  small  error  as  above  and 


~  2^P  )  <<  ♦  (  )  is  small  as  before  so  uie  need 

y  3  ^  ^ 

Cf)  «  1. 


For  the  smallest  particles  where  the  error  is  greatest 
r  3 

r  ft#  c  ;  thus  (,x/t)  1  (Pruppacher  and  Klett,  table  6.3). 

3 

There  is  therefore  considerable  error  in  the  cloud 
base  radius  calculated  for  these  small  particles  if  they 
have  a  very  large  insoluble  fraction.  However,  this  will 
not  introduce  significant  inaccuracies  to  the  droplet 
spectrum  calculated  in  the  cloud  ,  since  these  small 
particles  have  small  critical  radii  and  when  they  reach 
these  radii  the  error  is  small,  as  previously  shown.  For 
other  particles  the  error  is  smell  throughout  growth. 
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TO  FIGURES 

1  Measured  size  distribution  on  15th  May  1979  for  the 

period  1635-1650.  (A)  Kelly  probe,  (B)  ASSP. 

2  Case  I;  Acoustic  sounder  trace.  Dark  regions  represent 
echoes  from  thermal  plumes. 


3a  Case  I.  Observed  (A)  and  calculated  (B)  liquid  water 
contents  L. 

3b  Drop  concentration  N  and  percentage  P  of  histograms 
outside  N  ^  2/^ during  the  period  covered  by  case  I. 

U  Case  I.  A,  observed  size  distribution  over  the  period 
1331  to  13A6  aST  using  the  KP.  N  =  265  cm 
L  =  0.28  gm~^.  B,  calculated  spectrum  N  =  380  cm”^, 

L  =  0.31  gm"^. 

5  Case  II;  Acoustic  sounder  trace.  The  intense  layer-echo 
at  around  50-1Q0m  represents  the  base  of  the  inversion 
or  entrainment  interface. 

6  Drop  concentration  N,  percentage  P  of  histograms 
outside  N  2/f^,  and  mean  radius  r"  during  the  period 
covered  by  case  II 

7  Observed  (A)  and  calculated  adiabatic  (B)  liquid 
luatex  contents  L  during  period  covered  by  case  II. 

8  Case  II.  A,  measured  size  distributions  for  the  period 

0830  to  08^5  3ST.  N  =  132  cm  B,  calculated 

adiabatic  spectrum  N  =  212  cm~^. 

9  Case  II.  Size  distributions.  A,  observed  high  count 
spectrum,  B,  observed  low  count  spectrum.  A,B, 
averaged  over  15  minute  period  0300-0815  BST.  C, 
calculated  on  the  turbulent  mixing  model  of  Baker  and 
Lattiari  (19-’.1)  for  parameter  values  given  in  the  text. 
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Figure  10  Case  III,  acoustic  sounder  trace 

Figure  11  Observed  (A)  and  calculated  (0)  liquid  water  contents 


during  the  period  covered  by  case  III. 

Figure  12  Drop  concentration  N,  percentage  P  oF  histograms 

outside  N  and  mean  radius  7  during  the  period 

covered  by  case  III. 

Figure  13  Case  III.  Measured  size  distributions  for  the  period 
1236-1251  0ST.  A,  high  count  spectrum.  0,  low  count 
spectrum.  C,  calculated  on  the  turbulent  mixing  model 
Baker  and  Latham  (1931)  for  parameter  values  mentioned 
in  the  text.  A,  L  =  Q.2Ag  m"^;  a,C,  L  =  Q.ILg  m”^. 

Figure  1L  'Case  III.  A,  measured  size  distributions  for  the 
period  1257-1306  BST.  ¥  =  253  cm“^,  'C  =  0.47  gm"^ 

0,  measured  size  distributions  for  the  period  1622- 
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Figure  5 
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III.  II\!FRA-RED  ATTENUATION  IN  LIMITED  VISIBILITY  CONDITIOIMS 
Introduction 

Simultaneous  measurements  of  the  attenuation  of  infra¬ 
red  radiation,  cloud  droplet  size  spectra  and  cloud  liquid 
water  content  have  recently  been  made  at  the  summit  of  Great 
Dun  Fell.  Preliminary  analysis  of  these  data  suggests  that 
this  attenuation  at  about  11(jm  may  be  predicted  from  measurements 
of  cloud  liquid  water  content  by  means  of  the  Chylek  relation¬ 
ship  to  a  reasonable  degree  of  accuracy.  Good  temporal 
correlation  of  extinction  coefficient  values  derived  from 
transmissometer  measurements  with  those  derived  from  cloud 
liquid  water  contents  clearly  exists  and  further  analysis  is 
being  conducted  in  an  endeavour  to  establish  the  likely  errors 
in  the  extinction  coefficient  measurements  arising  from  factors 
such  as  multiple  scattering  within  the  transmission  path. 


Field  Measurements 


k6 


A  Barnes  Model  14  Infra-red  transmissometer  system 

available  under  a  collaborative  research  agreement  with  Royal 

Signals  and  Radar  Establishment  at  Malvern  was  operated  utilizing 

a  total  transmission  path  length  of  41.5  metres.  The  transmission 

measurements  were  made  utilizing  a  spike  filter  centred  at 

10.9am  0.3.im  though  the  instrument  is  also  capable  of  operating 

over  3-5um  or  3-13iam  radiation  bands.  It  uas  considered  that 

the  nominal  angular  field  of  view  of  3  milliradians  for  the 

radiation  receiver  uith  the  associated  very  small  detector 

size  rendered  insignificant  any  forward  scattering  corrections 

which  are  neglected  for  the  purposes  of  data  analysis. 

Simultaneous  measurements  of  cloud  droplet  spectra  were  made 

with  a  PMS  FSSP  instrument  which  is  capable  of  sizing  droplets 

of  3  to  47ijim  in  diameter.  In  addition,  direct  measurements  of 

cloud  liquid  water  content  were  made  with  an  optical  scattering 

instrument  developed  at  UMIST  which  possesses  a  sample  volume  of 

a  few  millilitres  and  a  path  length  of  approximately  0.2  metres. 

For  ground-based  usage,  the  FSSP  is  fitted  with  an  intake  horn 

and  an  aspirator  fan  which  provides  an  aspiration  velocity  of 

about  26  metres  per  second  through  the  sampling  region  and  an 

_  -1 

overall  volume  sampling  rate  of  8.14  ml  s  .  It  became  apparent 

_  'I 

during  these  field  studies  that  winds  in  excess  of  about  15  m  s 
exerted  a  ’ram*  effect  upon  the  intake  horn  of  the  FSSP  resulting 
in  significant  increases  in  the  volume  flow  rate  of  air  through 
the  instrument  and  consequent  overestimation  of  droplet  number 
concentrations.  Limited  investigation  of  this  phsnomenon 

_  1 

suggested  that  moderate  wind  speeds,  of  less  than  about  lOm  s  , 
did  not  significantly  influence  droplet  number  concentration 

measurements  but  that  overestimations  by  as  much  as  a  factor 

_ 

of  4  might  occur  with  wind  velocities  approaching  20m  s  .  The 


f 


optical  scattering  instrument  proved  suitable  for  ground-based 
operation  in  situations  uhere  cloud  liquid  water  contents  in 
excess  of  about  Q.05c  m“^  were  encountered.  Its  bandwidth 
of  about  IKHz  and  its  small  sample  volume  made  it  extremely 
responsive  to  small  scale  fluctuations  in  liquid  water  content. 


Preliminary  Data  Analysis 


Data  from  the  Barnes  transmissometer  was  recorded  in 
analogue  form  on  a  chart  recorder  and  values  of  extinction 
coefficient  for  a  given  cloud  sample  uere  calculated  by  taking 
20s  averages  of  the  chart  record  in  order  to  correspond  uith 
droplet  size  and  number  concentration  spectra  obtained  from  the 
FSSP  over  the  same  time  intervals.  These  spectra  uere  inte¬ 
grated  to  provide  a  aeries  of  measurements  of  cloud  liquid 
water  content.  Chylek  (1973)  suggested  that  extinction 
coefficient,  a-^,  could  be  related  to  liquid  water  content,  bl, 
by  means  of  the  relationship 
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''e  -  2pX 

where  C  is  a  constant  for  a  given  experimental  situation, 
p  is  the  water  density  and  X  is  the  wavelength  of  the  radiation. 
Thus,  values  of  extinction  coefficient  obtained  by  utilizing 
this  formula  could  be  compared  with  values  derived  from  direct 
transmission  measurements,  A  typical  sample  of  data  obtained 
by  this  method  is  plotted  in  Figure  1  and  exhibits  a  strong 
temporal  correlation  between  the  extinction  coefficient  curves 
provided  allowance  is  made  for  the  FSSP  sampling  errors.  Clearly 
the  FSSP  correction  factor  is  not  sufficiently  well  defined 
to  be  very  meaningful  here  and  demonstrates  the  need  for  more 
detailed  calibration  studies.  The  close  temporal  correlation 
is  nevertheless  remarkable  when  the  extended  sample  volume  of 
the  transmissometer  is  compared  with  that  of  the  FSSP,  and 
server  to  support  the  validity  of  this  experimental  approach. 

Figure  2  shows  a  comparison  of  extinction  coefficients 


from  attenuation  measurements  compared  with  ones  derived  from 
liquid  water  content  measurements  made  with  the  optical  scattering 
device  again  using  the  Chylek  formula.  Unfortunately,  the 
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period  analysed  does  not  exhibit  a  very  uide  range  of  values 
of  extinction  coefficient  but  nevertheless  a  reasonable  degree 
of  correlation  betueen  the  tuo  sets  of  measurements  may  be 
noted. 

Multiple  scattering  errors  arise  when  radiation  scattered 
out  of  the  transmitted  beam  is  re-scattered  into  the  beam  by 
interaction  with  another  particle.  This  multiple  scattering 
is  increased  as  the  number  of  scatterers  is  increased  whether 
by  a  greater  concentration  of  particles  per  unit  volume  or  by 
an  extended  path  length  and  its  effect  must  be  to  enhance  the 
transmission  of  radiation  through  the  medium  compared  with  that 
predicted  from  the  liquid  water  content  measurements.  More 
detailed  study  of  the  data  from  the  field  measurements  does 
indeed  support  this  hypothesis  and  a  simple  scheme  whereby  a 

I 

corrected  extinction  coef  f  icient  ,cjg ,  is  derived  from  the 
measured  coefficient  utilizing  a  formula  of  the  type 

Og  = 

where  l\l  is  the  droplet  number  concentration  and  a  is  a  constant 
for  a  given  experimental  configuration,  does  appear  to  yield 
much  closer  agreement  between  the  calculated  and  measured  values. 
A  sample  of  data  treated  in  this  manner  is  reproduced  in 
Figure  3.  This  approach  is  rather  rudimentary  at  present  and 
a  much  more  detailed  treatment  is  required. 


ERhoh  bah  for  B.  T 

8«nes  Polo, 3  / 


